Measurement of the Inelastic Proton-Proton Cross-Section at y/s = 7 TeV with the 

ATLAS Detector 



A first measurement of the inelastic cross-section is presented for proton-proton collisions at a 
center of mass energy \fs = 7 TeV using the ATLAS detector at the Large Hadron Collider. In a 
dataset corresponding to an integrated luminosity of 20 A 40 " 1 , events are selected by requiring hits 
on scintillation counters mounted in the forward region of the detector. An inelastic cross-section 
of 60.3 ±2.1 mb is measured for £ > 5 x 10 -6 , where £ = M x /s is calculated from the invariant 
mass, Mx, of hadrons selected using the largest rapidity gap in the event. For diffractive events 
this corresponds to requiring at least one of the dissociation masses to be larger than 15.7 GeV. 



PACS numbers: 13.85.Hd, 13.85.Lg, 12.40.Nn 
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Since the earliest days of particle physics, measure- 
ments of the total pp and pp cross-sections and their the- 
oretical understanding have been topics of much inter- 
est [l|. The cross-sections can not yet be calculated by 
quantum chromodynamics (QCD), and many approaches 
have been used to describe the existing measurements. 
General arguments based on unitarity, analyticity, and 
factorisation imply a bound (the Froissart bound 0, Q ) 
on the high-energy behaviour of total hadronic cross- 
sections. This bound is independent of the details of 
the strong interaction dynamics and states that the total 
cross-section can not rise faster than ln 2 (s), where y/s is 
the centre-of-mass energy. Recently it has been extended 
to the inelastic cross-section Existing experimental 
data [l[ show a rise in the hadronic cross-sections with s, 
but it is unclear whether the asymptotic behaviour has 
already been reached. With the data presented in this 
letter we shed further light on the high energy behaviour 
of the inelastic cross-section. 

The most common models that describe the data up 
to y/s — 1.8 TeV predict a rise of the cross-section 
with a simple power law (s Q (°) _1 where a(0) denotes 
the Pomeron-trajectory intercept) [H-@] or as a loga- 
rithm |8hfll|. Others employ QCD for aspects of the 
calculation [12| - [T3 |. However, while the phenomenolog- 
ical description of the existing data is largely adequate, 
there are significant uncertainties on the extrapolation to 
higher energies, partly due to a long-standing 2.7a dis- 
crepancy between the two highest energy collider mea- 
surements of the total pp cross-section by CDF [l5[ and 
by E811 (lfi. 

In this letter a measurement of the inelastic pp cross- 
section is presented using data taken by the ATLAS ex- 
periment [TtJ at the Large Hadron Collider (LHC) [l8[ at 
y/s = 7 TeV. The data considered were collected during a 
single 8-hour fill beginning March 31st, 2010, correspond- 
ing to an integrated luminosity of 20.3 ± 0.7 /ib -1 and 
a peak instantaneous luminosity of 1.2 x 10 27 cm _2 s _1 . 
The mean number of interactions per crossing in this 
fill is approximately 0.01. The analysis uses highly 
efficient scintillation counters to detect inelastic colli- 
sions. They are insensitive to diffractive dissociation 
processes in which the dissociation systems have small 
invariant masses, Mx- Their acceptance corresponds ap- 
proximately to £ = M x /s > 5 x 10~ 6 , equivalent to 



M x > 15.7 GeV for y/s = 7 TeV. The cross-section mea- 
surement presented here is restricted to this kinematic 
range. However, in order to compare the data with pre- 
vious measurements, an extrapolation of the cross-section 
is performed to the full £ range, £ > m 2 /s where m p is 
the proton mass. 

The ATLAS detector is described in detail else- 
where fl7| . The beam-line is surrounded by a tracking 
detector that uses silicon pixel, silicon strip, and straw 
tube technologies and is embedded in a 2 T magnetic 
field. The tracking system covers the pseudorapidity (l9j 
range \rj\ < 2.5. It is surrounded by electromagnetic 
and hadronic calorimeters covering \rj\ < 3.2, which are 
complemented by a forward hadronic calorimeter cov- 
ering 3.1 < \t}\ < 4.9. Minimum Bias Trigger Scintil- 
lator (MBTS) detectors, the primary detectors used in 
this measurement, are mounted in front of the endcap 
calorimeters on both sides of the interaction point at 
z = ±3.56 m and cover the range 2.09 < \rf\ < 3.84. 
Each side consists of 16 independent counters divided 
into two rings; the inner 8 counters cover the rapidity 
range 2.83 < |ry| < 3.84 and the outer 8 counters cover 
the range 2.09 < \r)\ < 2.83. Each individual counter 
spans 45° of the azimuthal angle (0), and 31 out of 32 
counters were operational. The luminosity is measured 
using a Cherenkov light detector, LUCID, which is lo- 
cated at z = ±17 m. The luminosity calibration has 
been determined during dedicated van der Meer beam 
scans to a precision of 3.4% [20l. l2lj] . 

Monte Carlo (MC) simulations are used to determine 
the acceptance of the event selection and to assess sys- 
tematic uncertainties. The detector response to the 
generated events is simulated using the ATLAS simula- 
tion [22| based on Geant4 [23[ , and both the simulated 
and data events are reconstructed and analysed with the 
same software. The Pythia6 [24], Pythia8 [25| and 
Phojet (2|| generators are used to predict properties of 
inelastic collisions. These generators distinguish between 
different processes that contribute to inelastic pp inter- 
actions: single dissociative (SD) processes, pp — > pX, in 
which one proton dissociates; double dissociative (DD) 
processes, pp — » XY, in which both protons dissociate 
with no net color flow between the systems X and Y; 
and non- diffractive (ND) processes in which color flow 
is present between the two initial-state protons. The 
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model by Schuler and Sjostrand (27j . used by Pythia6 
and Pythia8, predicts cross-sections of 48.5 mb, 13.7 mb 
and 9.3 mb for the ND, SD and DD processes, respec- 
tively. While the cross-sections used by Pythia6 and 
Pythia8 are identical, they differ in the modelling of the 
hadronic final state. Phojet predicts the correspond- 
ing cross-sections as 61.6 mb (ND), 10.7 mb (SD) and 
3.9 mb (DD). Due to differences in implementation of 
the interface between large £ diffractive (SD and DD) 
processes and ND processes in Pythia and Phojet the 
fractional contribution of these processes is a model de- 
pendent quantity. Phojet also includes a 1.1 mb contri- 
bution from central diffraction (CD), pp — > ppX, a pro- 
cess not implemented in Pythia, wherein neither proton 
dissociates but the Pomeron-trajectory exchange results 
in energy loss for the protons and the production of a 
central system of particles. The MC generators define 
the inelastic cross-section as the sum of these contribu- 
tions, and thus Schuler and Sjostrand (Phojet) predicts 
an inelastic cross-section of 71.5 mb (77.3 mb). Other re- 
cent predictions for this cross-section at y/s = 7 TeV are 
69 mb 0, 65-67 mb [13, 68 mb [3 and 60-75 mb [l|. 

The variable £ is defined at the particle level by divid- 
ing the final state particles into two systems, X and Y. 
The mean rj of the two particles separated by the largest 
pseudorapidity gap in the event is used to assign all par- 
ticles with greater pseudorapidity to one system and all 
particles with smaller pseudorapidity to the other [28| . 
The mass, Mx,y, of each system is calculated and the 
higher mass system is defined as X while the lower mass 
system is defined as Y. The variable £ is then given 
by £ = M x /s and it is bounded by the elastic limit of 
£ > rrip/s. Due to the limited MBTS detector accep- 
tance in ry, this measurement is restricted to the range 
£ > 5 x 10 -6 ; there is no restriction on My- Several mod- 
els are used for the dependence of the diffractive cross- 
sections on £. The Schuler and Sjostrand model has a 
relatively flat dependence on £, while the Phojet model 
predicts a slight decrease with decreasing £. Pythia8 
has several additional predictions for the ^-dependence of 
the diffractive cross-sections which are considered. Bruni 
and Ingelman [2!| predict a flat ^-dependence while Don- 
nachie and Landshoff (DL) [3(1, Berger et al. [3l|, and 
Streng [32| predict 
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where e = a(0) — 1. Values of e between 0.06 and 0.10, 
and of a' between 0.10 and 0.40 GeV -2 are considered for 
the DL model, a' is the slope of the Pomeron trajectory 
which is assumed to be linear such that a(t) — a(0) + a't. 
The DL model with e = 0.085 and a' = 0.25 GeV" 2 with 
Pythia8 fragmentation is the default model in this anal- 
ysis and the other models are used to assess uncertainties 
in the modelling of diffractive events. 

Experimentally the cross-section is calculated using 



where N is the number of selected events, Nbg is the 
number of background events, /^<5 X io- 6 i s the fraction of 
events that pass the event selection but have £ < 5 x 10~ 6 , 
J Ldt is the integrated luminosity, and ttrig and e se i are 
the trigger and offline event selection efficiencies in the se- 
lected £-range. Note that in this analysis £ is defined only 
at the particle level; the cut value at 5 x 10~ 6 was cho- 
sen such that the efficiency of the MBTS requirement is 
greater than 50% for any ^-values greater than 5 x 10~ 6 . 
In this measurement Nbg and e trig are determined di- 
rectly from the data. The MBTS individual counter ef- 
ficiencies in the MC simulation are tuned to match the 
observed efficiencies in data. Then e se i and /^<5 X io- 6 
are taken from the tuned MC simulation. In order to re- 
duce the uncertainties in the factors taken from MC sim- 
ulation, the relative diffractive dissociation cross-section, 
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Each of these steps is described in detail below. 

The MBTS functions as a trigger by determining the 
number of scintillation counters with a signal passing a 
leading-edge discriminator; in this analysis at least one 
trigger signal must be present. In the offline reconstruc- 
tion, the MBTS signals are fit to obtain the total charge 
and timing of the signal. The offline event selection re- 
quires at least two counters with a charge larger than 

0. 15 pC. This threshold is set to be well above the noise 
level, which is well described by a Gaussian centred at 
zero of width 0.02 pC. This inclusive sample contains 
1,220,743 data events. In order to constrain the diffrac- 
tive components a subset of events is selected, the single- 
sided sample, which contains events that have at least 
two hits on one side of the MBTS detector and no hits 
on the opposing side (in z). In the data 122,490 single- 
sided events are observed. 

Backgrounds arise from beam-related interactions, 
such as collisions of the beam with gas particles in the 
beam-pipe or with material upstream from the detector, 
and slowly-decaying, collision-induced radiation termed 
"afterglow" [2l|. Additionally, instrumental noise and 
cosmic rays provide backgrounds which were studied and 
found to be negligible for this analysis. The beam-related 
backgrounds are determined using the number of selected 
events collected in this fill with the non-colliding bunches, 

1. e. when only one proton bunch was passing through AT- 
LAS [33j], normalised by the ratio of the number of pro- 
tons in the colliding to the non-colliding bunches. The 
single-sided selection contains 422±28 background events 
and the inclusive sample contains Nbg = 1,574 ± 54 
background events, corresponding to 0.3% and 0.1% of 
the total samples, respectively. In addition, there is an 
in-time afterglow component due to the scattering of sec- 
ondary low-energy particles produced in the same colli- 
sion event which can give additional hits, causing low- 
activity events to migrate into the selected event sample. 
This contribution is evaluated to be at most 0.4% for the 
inclusive, and 3.6% for the single-sided, samples by ex- 
amining the asymmetry of the absolute timing measure- 
ment of the MBTS counters. We conservatively assume 
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a 100% uncertainty on both background sources which 
covers any residual impact of the afterglow on the back- 
ground subtraction, any uncertainty in the beam current 
measurements and the uncertainty due to in-time after- 
glow. The resulting overall uncertainty on the number of 
background events Nbg is given by the quadratic sum of 
the two components and is 0.4%. 

The trigger efficiency of the MBTS detector with re- 
spect to the offline requirement, e tr i g , is measured to 
be 99.98t°;52% (statistical errors) using events triggered 
randomly on colliding beams. The systematic uncer- 
tainty on etrig is determined using a second, independent 
trigger as reference. The difference between the two effi- 
ciency determinations leads to a 0.1% uncertainty on the 
cross-section measurement. 

The data and MC simulation agreement in the MBTS 
counter response is checked using other detector subsys- 
tems with overlapping r\ ranges: charged particles recon- 
structed by the tracking detector (2.09 < \q\ < 2.5), and 
calorimeter showers in the inner wheel of the electromag- 
netic calorimeter (2.5 < \t)\ < 3.2) and in the forward 
calorimeter (3.1 < \r/\ < 3.84). The efficiency with re- 
spect to a track (calorimeter energy deposit) to have a 
signal above the 0.15 pC threshold in the outer (inner) 
counters is on average 98.5% (97.5%) for the data and a 
constant 99.4% (98.7%) in the MC simulation. The indi- 
vidual counter efficiencies deviate by up to 2.0% (2.5%) 
from the average in the data. The MC simulation is cor- 
rected to match the data efficiency and the maximum 
variations in the counter responses are considered as a 
systematic uncertainty. This results in a 0.1% uncer- 
tainty on the cross-section measurement. 

The offline selection efficiency, e se ;, depends on the 
amount of material traversed by particles before hitting 
the MBTS detector. The rate of photons (primarily from 
7T° decays) converting to electrons which are subsequently 
detected by the MBTS increases with additional mate- 
rial, resulting in an increase of e se ;. Second order ef- 
fects arise from charged particles scattering out of the 
MBTS acceptance region (decreasing e se ;), or charged 
particles scattering into the acceptance region (increas- 
ing /{<5xio- 6 )- Within the tracking volume (|?7| < 2.5) 
the material distribution has been studied using conver- 
sion electrons and K® — > ir + 7r~ decays, and is known 
to better than ±5% in the central region of the detec- 
tor and to ±30% for 2.2 < \r)\ < 2.5 33]. In the region 
\rj\ > 2.5 the material is dominated by the cooling and 
electrical services to the silicon pixel detector, and an un- 
certainty of ±40% is assumed. This is validated in-situ 
using the fraction of events where we observe significant 
energy in the forward calorimeters but no signal (above 
noise) in the MBTS detector. The resulting systematic 
uncertainty on the cross-section is 0.2%. 

Misalignments of the MBTS detector with respect to 
the nominal centre of the detector could change the event 
selection efficiency for a particular value of £. Misalign- 
ments of up to 10 mm were considered and found to have 
a negligible impact. A misalignment of 10 mm is conser- 



vative compared to the survey precision and any known 
misalignments within the ATLAS experiment [17] . 

The fractional contribution of diffractive events, fu, 
is constrained by the ratio of single-sided to inclusive 
events, R ss - The MC generators predict that less than 
1% of the ND process pass the single-sided event selec- 
tion, while 27 — 41% of the SD and DD processes pass 
the single-sided selection. For all models the inclusive 
sample is dominated by ND events, therefore the ratio of 
single-sided to inclusive events is sensitive to the relative 
fraction of diffractive events. 

The measured R ss in the data is R ss — [10.02 ± 
0.03(stat.) toi( s y s t-)]%; where the systematic error in- 
cludes the uncertainties on the backgrounds, the MBTS 
response and the material. 




Figure 1: The ratio of the single-sided to inclusive event sam- 
ple Rss as a function of the fractional contribution of diffrac- 
tive events to the inelastic cross-section fr>. The data value 
for Rss is shown as the horizontal line with its systematic un- 
certainties (grey band). Also shown are predictions of several 
models as a function of an assumed value of /d- The default 
fr> value (32.2% for all models but Phojet which is 20.2%) 
is indicated by the markers. 



Fig.rjjcompares the observed value of R ss to the predic- 
tions of several models as a function of fn- The intersec- 
tion of the R ss value measured in data with the prediction 
is used as the central value of fu for each model. The 
systematic uncertainty on fjj is determined by the maxi- 
mum and minimum values consistent with the la uncer- 
tainty on the data when varying the double- to single- 
dissociation event ratio between and 1. The resulting 
value using the default DL model is fo = 26.9ti'g%- 

The acceptance calculation relies on the MC genera- 
tors to provide an adequate description of the particle 
multiplicity in the acceptance region. The validity of the 
MC description is assessed by examining the hit multi- 
plicity in the MBTS detector in the inclusive and single- 
sided event samples as shown in Fig. [2] While none of 
the generators gives a perfect description, the data lie 
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between the models at low multiplicity which is most 
important for the measurement. The default DL model 
describes the single-sided sample well, giving confidence 
in the diffractive modelling. We use the difference in 
the MC correction factor determined with Pythia8 and 
Pythia6 as the uncertainty due to the fragmentation 
model, leading to a 0.4% uncertainty on the cross-section. 
The maximum difference between the default DL model 
and all other models is taken as the uncertainty due to 
the underlying £ distribution. Variations of a' have a 
negligible effect on the acceptance. Among all the models 
considered, the Phojet model gives the largest difference 
in the correction factor, leading to a 0.4% uncertainty on 
the cross-section. 



Source 


Uncertainty (%) 


Trigger Efficiency 


0.1 


MBTS Response 


0.1 


Beam Background 


0.4 


fo 


0.3 


MC Multiplicity 


0.4 


^-Distribution 


0.4 


Material 


0.2 


Luminosity 


3.4 


Total 


3.5 



Table I: Sources of systematic uncertainty and their effect on 
the cross-section measurement. 
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Figure 2: The MBTS multiplicity distribution in the data 
(fill ed ci rcles) compared with MC expectations for the inclu- 
sive (a) and single-sided (b) samples for several MC models 
(histograms) using the fitted fo values. The band around the 
data indicates the systematic uncertainty due to the MBTS 
detector response and the amount of material in front of the 
MBTS detector. 



The final result for the measured inelastic cross-section 
is calculated using the default DL model of e = 0.085 
and a' = 0.25, which yields f D = 26.9%, e se i = 98.77%, 
and f (<5xW -e = 0.96%. Together with e trig = 99.98%, 
N = 1,220,743, N BG = 1,574 and / Ldt = 20.25 /ib" 1 
this results in cr meZ (£ > 5 x 10" 6 ) = 60.3 ± 0.05(stat) ± 
0.5(syst) ± 2.1(lumi) mb. The systematic uncertainty in- 
cludes all contributions discussed above and listed in Ta- 
ble |U the dominant uncertainty arises from the luminos- 
ity calibration and is quoted separately. 

The measurement is compared to the predictions 
in Figure [3] and Table HU The predictions by the 
Schuler-Sjostrand model (66.4 mb) and the Phojet 
model (74.2 mb) are both higher than the data. The 
prediction of 51.8-56.2 mb by Ryskin et al. [HJ is slightly 
lower than the data. 

In order to compare with previous measurements and 
analytic models, the fractional contribution to the in- 
elastic cross-section of events passing the £ > 5 x 10~ 6 
cut is determined from the models and used to extrap- 
olate the measurement to the full inelastic cross-section. 
This fraction is 87.3% for the default model of DL with 
e = 0.085 and a' = 0.25. The other models considered 
give fractions ranging from 96%(Phojet) to 86%(DL 
with e = 0.10). Recent calculations also yield values 
between 79% and 84% [Hj]. Thus 87% is taken as the 
default value for this fraction and an uncertainty of 10% 
is taken due to the extrapolation uncertainty on the in- 
dependence. The resulting inelastic cross-section value 
is crenel = 69.4 ± 2.4(exp.) ± 6.9(extr.) mb where exp in- 
cludes the statistical and experimental systematic errors, 
including the luminosity uncertainty. 

This result is shown in Fig. 3 compared to several theo- 
retical predictions and a variety of data at lower yfs. The 
measurement within the kinematic range £ > 5 x 10~ 6 
is significantly lower than the predictions of Schuler and 
Sjostrand and Phojet. The extrapolated value agrees 
within the large extrapolation uncertainty with the pre- 
dictions from Pythia, which uses a power law depen- 
dence on y/s. It also agrees with Block and Halzen [l(| 
(which has a logarithmic yfs dependence) , and with other 
recent theoretical predictions that vary between 60 and 
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Figure 3: The inelastic cross-section versus y/s. The ATLAS 
measurement for £ > 5 x 10 -6 is shown as the red filled circle 
and compared with the predictions of Schuler and Sjostrand 
and Phojet for the same phase space. Data (filled circles 
for pp data and unfilled circles for pp data) from several ex- 
periments are compared with the predictions of the pp inelas- 
tic cross-section from Schuler and Sjostrand [27j ] (as use d by 
Pythia), by Block and Halzen [l(| and by Achilli et al. [Tjf. 
An extrapolation from the measured range of £ > 5 x 10~ 6 to 
the full inelastic cross-section using the acceptance of 87±10% 
is also shown (blue filled triangle). The experimental uncer- 
tainty is indicated by the error bar while the total (includ- 
ing the extrapolation uncertainty) is represented by the blue 
shaded area. 



72 mb [12h14| . It should be stressed that this extrapo- 
lation relies on the prediction of the ^-dependence of the 
cross-section. 

The measurement and a variety of theoretical predic- 
tions are also summarised in Table [TT1 

In conclusion, a first measurement of the inelastic 
cross-section has been presented for pp collisions at \fs = 
7 TeV with a precision of 3.5%. The measurement is lim- 
ited to the kinematic range corresponding to the detector 
acceptance: £ > 5 x 1CP 6 . Phenomenological predictions 
for both a power law dependence and a logarithmic rise 
of the cross-section with energy are consistent with the 
measurement. 
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5 x 10~ 6 , and for the full kinematic range, £ > m 2 p /s. The 
experimental uncertainty (exp.) includes the statistical, sys- 
tematic and luminosity uncertainties. The extrapolation un- 
certainty (extr.) only applies to the full kinematic range and 
is listed separately. 



ICYT, Chile; CAS, MOST and NSFC, China; COL- 
CIENCIAS, Colombia; MSMT CR, MPO CR and VSC 
CR, Czech Republic; DNRF, DNSRC and Lundbeck 
Foundation, Denmark; ARTEMIS, European Union; 
IN2P3-CNRS, CEA-DSM/IRFU, France; GNAS, Geor- 
gia; BMBF, DFG, HGF, MPG and AvH Foundation, 
Germany; GSRT, Greece; ISF, MINERVA, GIF, DIP and 
Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, 
Japan; CNRST, Morocco; FOM and NWO, Netherlands; 
RCN, Norway; MNiSW, Poland; GRICES and FCT, 
Portugal; MERYS (MECTS), Romania; MES of Rus- 
sia and ROSATOM, Russian Federation; JINR; MSTD, 
Serbia; MSSR, Slovakia; ARRS and MVZT, Slovenia; 
DST/NRF, South Africa; MICINN, Spain; SRC and 
Wallenberg Foundation, Sweden; SER, SNSF and Can- 
tons of Bern and Geneva, Switzerland; NSC, Taiwan; 
TAEK, Turkey; STFC, the Royal Society and Lever- 
hulme Trust, United Kingdom; DOE and NSF, United 
States of America. 

The crucial computing support from all WLCG part- 
ners is acknowledged gratefully, in particular from 
CERN and the ATLAS Tier-1 facilities at TRIUMF 
(Canada), NDGF (Denmark, Norway, Sweden), CC- 
IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF 
(Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Tai- 
wan), RAL (UK) and BNL (USA) and in the Ticr-2 fa- 
cilities worldwide. 



6 



[1] K. Nakamura et al. (Particle Data Group), Chapter 41, [19] 
pg.385-401, J. Phys. G 37 (2010) 075021. and references 
therein. 

M. Froissart, Phys. Rev. 3 123 (1961) 1053-1057. 
A. Martin, Nuovo cimcnto 42 (1966) 930. 
A. Martin, Phys. Rev.D 80 (2009) 065013. 
P.D.B. Collins, An introduction to Regge theory and high 
energy physics, Cambridge University Press 1977, and 
references therein. 

A. Donnachie and P.V. Landshoff, Phys. Lett. B 296 [20 
(1992). [21 
A. Donnachie and P.V. Landshoff, Phys. Lett. B 595 
(2004) 393. 

M.M. Block, K. Kang and A.R. White, Int. J. Mod. Phys [22 
A7 (1992) 4449, and references therein. [23 
M. M. Block and R. N. Cahn, Rev. Mod. Phys. 57 (1985) 
563; M. M. Block, Physics Reports 436 (2006) 71-215. 
M. M. Block and F. Halzen, arXiv:1102.3163. [24 
J. R. Cudell, K. Kang, S. K. Kim, Phys. Lett. B 395 
(1997) 311-331. [25 
M. Ryskin, A. Martin, V. Khoze, arXiv:1102.2844, sub- 
mitted to Eur. Phys. J. C; private communication. [26 

E. Gotsman, E. Levin, U. Maor, J. S. Miller, Eur. Phys. 
J. C57 (2008) 689-709; E. Gotsman, E. Levin, U. Maor, [27 
arXiv:1010.5323. 

A. Achilli et al, CERN-PH-TH/2011-018, [28 
arXiv:1102.1949, submitted to JHEP. 

F. Abe et al. (CDF Collaboration), Phys. Rev. D 50 [29 
(1994) 5550-5561. [30 
C. Avila et al. (E811 Collaboration), Phys. Lett. B 445 
(1999) 419-422. [31 
ATLAS Collaboration, JINST 3 (2008) S08003. [32 
L. Evans, (ed.), P Bryant, (ed.), LHC Machine, JINST [33 
3 (2008) S08001. 



ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre 
of the detector and the 2-axis coinciding with the axis 
of the beam pipe. The a>axis points from the IP to the 
centre of the LHC ring, and the y axis points upward. 
Cylindrical coordinates (r, <f>) are used in the transverse 
plane, cf> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 
6 as r) = -lntan(0/2). 

S. van der Meer, CERN-ISR-PO-68-31, 1968. 
ATLAS Collaboration, ATLAS-CONF-2011-011; AT- 
LAS Collaboration, arXiv:1101.2185vl, accepted by 
EPJC. 

ATLAS Collaboration, Eur. Phys. J. C70 (2010) 823-874. 
S. Agostinelli, et al., GEANT4 Collaboration, GEANT4: 
A simulation toolkit, Nucl. Instr. Meth. A 506 (2003) 
250303. 

PYTHIA version 6.421; T. Sjostrand, S. Mrenna, P. 
Skands, JHEP 05 (2006) 026. 

PYTHIA version 8.135; T. Sjostrand, S. Mrenna, P. 

Skands, Comput. Phys. Commun. 178 (2008) 852-867. 

R. Engel, Z. Phys. C 66 (1995) 203-214; R. Engel, J. 

Ranft, Phys. Rev. D 54 (1996) 4244-4262. 

G. Schuler and T. Sjostrand, Phys. Rev. D 49 (1994) 

2257-2267. 

C. Adloff et al. (HI Collaboration), Z. Phys. C 74 (1997) 
221. 

P. Bruni and G. Ingelman, Phys. Lett. B311 (1993) 317. 
A. Donnachie and P. V. Landshoff, Nucl. Phys. B244 
(1984) 322. 

E. L. Berger et al., Nucl. Phys. B286 (1987) 704. 

K. Streng, CERN-TH.4949 (1988) 

ATLAS Collaboration, Phys. Lett. B 668 (2010) 1. 



